OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Molecular Signatures of the Evolving Immune Response 

in Mice following a Bordetella pertussis Infection cros^k 

Rene H. M. Raeven^'^% Jolanda Brummelman^, Jeroen L. A. Pennings'*, Olaf E. M. Nijst\ Betsy Kuipers^, 
Laura E. R. Blok\ Kina Helm^ Elly van Riet\ Wim Jiskoot^ Cecile A. C. M. van Els^ Wanda G. H. Han^ 
Gideon F. A. Kersten^'^ Bernard Metz^ 

1 Intravacc, Brlthoven, The Netherlands, 2 Division of Drug Delivery Technology, Leiden Academic Centre for Drug Research, Leiden, The Netherlands, 3 Centre for 
Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands, 4 Centre for Health Protection (GZB), National 
Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands 



Abstract 

Worldwide resurgence of pertussis necessitates the need for improvement of pertussis vaccines and vaccination strategies. 
Since natural Infections Induce a longer-lasting Immunity than vaccinations, detailed knowledge of the Immune responses 
following natural Infection can provide important clues for such Improvement. The purpose was to elucidate the kinetics of 
the protective immune response evolving after experimental Bordetella pertussis (6. pertussis) infection in mice. Data were 
collected from (I) Individual analyses, i.e. microarray, flow cytometry, multiplex Immunoassays, and bacterial clearance; (li) 
twelve time points during the Infection; and (lil) different tissues involved In the Immune responses, i.e. lungs, spleen and 
blood. Combined data revealed detailed Insight in molecular and cellular sequence of events connecting different phases 
(Innate, bridging and adaptive) of the Immune response following the infection. We detected a prolonged acute phase 
response, broad pathogen recognition, and early gene signatures of subsequent T-cell recruitment in the lungs. Activation 
of particular transcription factors and specific cell markers provided insight into the time course of the transition from innate 
towards adaptive immune responses, which resulted In a broad spectrum of systemic antibody subclasses and splenic Thl/ 
Thl7 memory cells against 6. pertussis. In addition, signatures preceding the local generation of Thl and Thl 7 cells as well as 
IgA In the lungs, considered key elements in protection against B. pertussis, were established. In conclusion, molecular and 
cellular Immunological processes in response to live 6. pertussis infection were unraveled, which may provide guidance in 
selecting new vaccine candidates that should evoke local and prolonged protective immune responses. 
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Introduction 

The gram-negative bacterium Bordetella pertussis, the causative 
agent of whooping cough (pertussis), accounted for high mortality 
rates among infants in the pre-vaccrne era. Whole cell and 
acellular vaccines have drastically decreased the number of cases 
[1]. However, recently resurgence of pertussis in the vaccinated 
population has been observed in the USA [2,3] and in European 
countries [4,5]. Currently pertussis is stUl endemic, ranking in the 
top-ten of vaccine preventable diseases worldwide, according to 
the WHO. Next to improved diagnostics and public awareness, 
also strain adaptation and waning immunity are thought to be 
responsible for this increase in disease cases [6]. Thus, a call for 
vaccines with improved efficacy is evident. 

As a main feature of the innate immune response upon B. 
pertussis infection in mice the recognition of lipopolysaccharide 
(LPS) by TLR4 has been acknowledged [7,8]. This activation of 
TLR4 by B. pertussis leads to up-regulation of cytokine gene 
expression and recruitment of neutrophils into the lungs [9,10]. It 



was found in animal models that B. pertussis infection results in 
formation of T-helper (Th) 1 and Thl7 cells [11—13]. Since 
immunity induced by natural infections provides faster clearance 
upon reinfection and is longer lasting compared to both acellular 
and whole cell pertussis vaccination [14,15], immune mechanisms 
induced upon infection or vaccination have been compared. In 
human and murine studies, immunization with whole cell or 
acellular pertussis vaccines results predominantiy in a Thl or a 
Th2 response, respectively [11,16]. In addition, in both the 
intramuscular (human) or subcutaneous (mice) administered 
acellular and whole cell pertussis vaccines, the humoral response 
is characterized by systemic IgG [17,18], while mucosal immune 
responses seem absent. Despite the absence of direct evidence for 
correlates of protection against B. pertussis, both Thl and Thl 7 
type 004""" T-ceUs as well as IgA-producing B-cells seem to play an 
important role in a protective immune response against B. 
pertussis [12,15,19-22]. Since these responses are not induced 
by the currently available vaccines [23], but more resemble 
responses induced by infection, a better understanding of the 
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immune response induced in the host following B. pertussis 
infection is needed. Despite knowledge on particular elements of 
the immune response generated by a B. pertussis infection, littie is 
known about the kinetics and sequential relation of these elements. 
For this, systems biology can be an important tool, as was shown 
for tuberculosis and influenza infection [24—26] . 

Here, systems biology was applied to elucidate molecular and 
ceUular events in the difierent phases of the immune response after 
primary B. pertussis infection in a murine model. To this end, 
innate and adaptive immune responses were investigated over a 
period of 66 days post infection. Gene expression profiles in spleen 
and lungs, cytokine profiles in sera, and cellular composition of the 
spleen were determined at twelve time points. Furthermore, 
ceUular and antibody mediated immune responses against B. 
pertussis were investigated. Herewith, we revealed a chronological 
cascade of immunological processes consisting of recognition, 
processing, presentation and clearance of B. pertussis. The 
extensive insights into the immune response upon B. pertussis 
infection generated in this study may serve as a solid base for 
future research on pertussis vaccines and vaccination strategies. 

Results 

Lung clearance of infected mice 

The presence of B. pertussis in lungs of mice was examined 
during a period of 28 days post infection (p.i.), providing the 
benchmark for this study (Figure lA). Therefore, mice were 
intranasaUy infected with B. pertussis using a dose of lO'' colony- 
forming units (cfu). Approximately of the intranasal dose was 
traceabk" in the lungs of mice 2 hours p.i. The number of bacteria 
remained fairly constant for one day, and increased from the 
second day to a maximum 7 days p.i. (lO' cfu). Subsequentiy, a 
decrease in the number of bacteria was observed and complete 
clearance in 2 out of 3 mice was achieved 28 days p.i. To 
determine whether single intranasal infection with B. pertussis 
leads to protection, mice were reinfected 56 days after primary 
infection (Figure IB). A similar number of viable bacteria was 
detected 4 hours p.i. in lungs of both reinfected and naive mice. 
Reinfected mice were able to clear B. pertussis from the lungs 
within 2 days p.i., whereas naive mice showed a similar pattern as 
observed before. In conclusion, naive mice can clear B. pertussis 
from the lungs in about 28 days. Furthermore, mice previously 
infected with B. pertussis had developed sterilizing immunity, 
which clears the lungs in two days. 

Gene expression in lung tissue 

The gene expression in lung tissue of infected mice was 
monitored over a period of 28 days. In total 558 genes of the 
genome were dififerentially regulated (jti-valueSO.OOl, fold ratio 
(FR) >1.5): 446 genes were up-r(^gulated and 112 genes were 
down-regulated (Figure 2). A compk'te list of gene names is given 
in Table SI. Gene expression profiles changed already 4 hours p.i. 
with largest difierences at 14 days p.i. (377 responsive genes). The 
number of genes with an altered expression had declined 21 and 
28 days p.i. However, the expression of many genes had not yet 
returned to naive basal levels. Principal component analysis 
confirmed the largest differences at 14 days p.i. and the decline at 
21 and 28 days p.i. (Figure 3). Interestingly, the time point with 
the least changes in gene expression (2 days p.i.) preceded the 
increase in cfu in the lungs (Figure 1). 

To determine which biological processes were involved, 
functional analysis of all 558 genes was performed using DAVID. 
Over-representation analysis (ORA) (Benjamini-corrected p-vul- 
ueS0.05) in GO-BP and KEGG databases resulted in respectively 



122 terms and 20 pathways significanfly enriched. A selection of 
these terms is presented in Figure 4. The selection is based on 
immunological relevance and by exclusion of overlapping terms. 
Acute phase response and chemotaxis were obserx^ed right after 
infection (4 hours p.i.) (Figure 2). Subsequentiy, innate immune 
response, phagocytosis, antigen processing and antigen presenta- 
tion were initiated (7 days p.i.). Involvement of B-ceU mediated 
immunity in this stage is mainly demonstrated by the expression of 
complement related genes. Eventually, the gene expression 
showed initiation of the adaptive immune response, such as 
activation of B-cells and T-cells, B-cell receptor (BCR) signaling, 
humoral immune response (complement activation) and immuno- 
globulin-mediated immunity from 10 until 28 days p.i. 

The gene expression data was compared with the BioGPS 
database to identify the influx, presence or activation of particular 
immunological cells in the lungs (Figure SI). Sixty-one genes, 
which are predominantly expressed in macrophages, suggest two 
events: (i) triggering of alveolar macrophages 4 hours p.i. and (ii) 
recruitment of macrophages in the lungs 7 days p.i. The influx of 
macrophages was observed on cellular level using fluorescence 
microscopy [27]. In addition, the increased expression of 29 genes 
was attributed to the presence of neutrophils in the lungs 4 days 
p.i. Moreover, data suggests altered expression of 32, 48 and 17 
genes of B-ceUs, dendritic ceUs (DCs) and mast cells, respectively, 7 
days p.i., and 19 T-cefl genes 14 days p.i. 

DifiFerentiaUy regulated genes in lung tissue after the infection 
with B. pertussis were classified according to function, pathway 
and cell type. For ten different groups details are des[:ribed below. 

Cytokines. Twenty-one cytokine-encoding genes were differ- 
entially expressed in the lungs during B. pertussis infection 
(Figure 5). Sixteen cytokines are associated with a chemotactic 
function. The other five genes {Flt3l, Ifng, III 7a, 1117 f, 1116) 
reflect cytokines involved in the activation of T-cells, such as Thl 
and Thl 7 cells. 

Membrane receptors. Forty-six membrane receptors were 
identified with altered gene expression (Figure 5). The receptors 
were divided in eight groups: cytokine receptors, C-type lectins 
(CLRs) [28], Fc-receptors, epithelial cell specific receptors, 
including the IgA transporter (pIgR), myeloid-associated immu- 
noglobulin-like receptors, semaphorins [29], lymphocyte receptors 
(e.g., Pdll) [30,31] and the mucosal-homing receptor (Itgh7). 
Other membrane receptors are collected in one group (Figure 5). 

Gene expression profiles of both cytokines and cytokine 
receptors enabled early prediction of immune cell recruitment 
towards the lungs on transcriptomic level in this study. Increased 
expression of Cxcll3 and Ccl2() was followed by the up-regulation 
of receptors Cxcr5 and Ccr6. These chemokines lead to the 
recruitment of CXCR5+ B-ceUs, CXCR5+ T-ceUs, CCR6+ DCs 
and CCR6-I- Thl 7 ceUs into the lungs [32-34]. Likewise, tiie 
expression of Cxd9, CxcllO and Cxclll, induced by interferon 
gamma (IFNy), is associated with r('cruitmciit of CXCR3+ T-cells 
or CXCR3+ pDCs [35]. 

Antunicrobial peptides (AMPs). AMPs have defensive 
properties against bacterial infections, recruit immune cells and 
model subsequent immune responses [36]. Gene expression of 
four AMPs was elevated for 21 days p.i. (Figure 6). These four 
AMPs recruit neutrophils or protect against respiratory pathogens 
[37-40]. Interestingly, the chemokine Cxcll7, primarily expressed 
in the lungs, was described as an antimicrobial mucosal chemokine 
that recruits DCs and CD 14"^ monocytes and protects against 
pathogenic bacteria [41,42]. 

Acute phase response. Primary reaction to an infection is 
the acute phase response, which leads to the release of pro- 
inflammatory cytokines and activation of inflammatory cells. 
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Figure 1. Lung clearance of naive and reinfected mice after B. pertussis infection. (A) Number of colony forming units (cfu) in challenge 
culture (C.C.) was confirmed before challenge. All other cfu were determined in lung of challenged mice (mean n = 3). A large fraction of the original 
infection dose was traceable in the lungs of mice 2 hours p.i. Bacteria were able to colonize and multiply approximately 100-fold at 7 days p.i. After 
one week, the mice were able to clear bacteria, which resulted in cleared lungs (2 out of 3 mice) at the last time point (day 28 p.i.). (B) Reinfection was 
performed at 56 days after primary infection and the number of cfu were counted after 4 hours p.i. While reinfected mice were able to clear 6. 
pertussis from the lungs within 2 days p.i., naive mice showed a similar pattern as observed in Figure 1A. 
doi:1 0.1 371 /journal.pone.01 04548.g001 



Strong up-regulation of 19 genes involved in the acute phase 
response was detected at 4 hours p.i. Remarkably, the gene 
expression of the acute phase response remained liighly expressed 
at least until day 28 p.i. (Figure 6). 

Pathogen recognition. Enhanced gene expression suggested 
activation of TLR4, TLR9, DAI and CLR sigiialiiig pathways, 
implying pathogen recognition by several ligands of B. perlussis 
(Figure 6). The induction of Md-1 and transcription factor IrfS 
genes suggested TLR4-signaling upon sensing of LPS [43]. Up- 
regulation of Ctsk and the downstream transcription factor Irf7 
assumes activation of TLR9 by sensing of unmethylated CpG 
sequences in bacterial DNA [44] . In addition, the cytosolic DNA 
sensor (DAI) was induced [45]. Elevated gene expression of several 
CLRs was noticed, such as Mgll and the mouse DC-SIGN-related 
proteins Signr? and SignrS [46]. Mgll recognizes A^-acetylga- 
lactosamine (GalNAc), which is present in B. pertussis. Mgll has 
recently been linked to B. pertussis mast cell interaction [47]. 

Transcription factors. Ten transcription factors were dif- 
ferentially expressed in the lungs (Figure 6). These include genes of 
factors Pu.l, Shp-1, NfiB, Syk, IrfS, Irfl and IrfS. All of these 
transcription factors have been designated as regulators of 
immunological processes, thereby defming the direction of the 
host's immune response upon bacterial infection [48-54]. 

Phagocytosis. Six genes representing subunits of NADPH 
oxidase (NOX) were activated between 4 hours and 28 days p.i. 
(Figure 6). The activation suggests the production of reactive 
oxygen species, which are involved in a phagocyte respiratory 
burst to assist in the killing of bacteria during phagocytosis [55]. 

Antigen presentation. Twenty-eight genes representing an- 
tigen presentation were up-regulated between 7 and 28 days p.i. 
(Figure 6). Results indicate that both antigen presentation via 
MHC class I and class II were included. This suggests antigen 
presentation to both CD4"'" T-helper cells and cytotoxic CD8^ T- 
ceUs (CTLs). In line with these findings, the up-regulation of Cd4 
and Clstu gene expression in lung tissue was observed, which may 
indicate presence of Th-cells [56] and CTLs [57], respectively 
(Figure SI). 



Complement activation. Twelve genes of the complement 
system were more abundandy expressed between 7 and 28 days 
p.i. (Figure 6). Increased gene expression was observed for the 
complement components C 1 and C4, which belong to the classical 
pathway. In contrast, reduced gene expression of the C 1 -inhibitor 
{Serpingl) was observed day 4 p.i. Furthermore, strong up- 
regulation of C/6 suggests activation of the alternative complement 
pathway. 

B-cell receptor (BCR) signaling. Fifteen genes of the BCR 
signaling pathway were up-regulated between 7 and 28 days p.i. 
(Figure 6). BCR activation was detected by elevated gene 
expression of eight downstream proteins [Blnh, Cardll, Grb2, 
Pik3cd, Plcg2, Prkcb, Rac2 and Vavl). Furthermore, the gene 
expression of the BCR signaling pathways indicates a balanced B- 
ceU response because of an interplay of stimulation by Cd79a and 
Ifitml, and inhibition by Fcgrllh, Lilrb3 {Pirb), Cd22, Cd72 and 
Shp-1. These co-inhibitors prevent overstimulation of B-ceUs [58]. 

Immunoglobulins. Gene expression associated with anti- 
body formation was observed after 4 days p.i. and was increased 
up to 28 days p.i. (Figure 6). The microarray data revealed gene 
expression of seven kappa (k) light chain genes, heavy chain Ig 
alpha chain C {Ihg-vj358), Ighm and Ighg, indicating the 
formation of IgA, IgM and IgG. 

In summary, approximately 40% of the genes expressed in the 
lungs could be annotated according to function. The chronology 
of immunological processes comprises broad pathogen recogni- 
tion, prediction of cell recruitment by cytokine and receptor 
expression, MHC class I and II antigen presentation, and markers 
for both antibody and cell-mediated immunity. 

Cytokine profiles in mouse sera 

The concentration of thirty-three cytokines was determined in 
mouse sera by ELISA and multiplex immunoassay (MIA). B. 
pertussis infection induced nine cytokines: CCLll, CCL20, 
CXCL9, CXCLIO, G-CSF, IL-la, IL-6, IL-9 and IL-17A 
(Figure 7). The pro-inflammatory cytokines, IL-la and IL-6, are 
both involved in the acute phase response. The concentration of 
IL-la was significantly increased 1 day p.i. Enhanced concentra- 
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Figure 2. Pulmonary gene expression profiles as function of time after B. pertussis infection. Fold changes In expression were calculated 
compared to naive mice and significant gene expression results (FRal.5, p-value<0.001) are visualized as heatmap (mean of n = 3). In total, 558 
genes were found to be differentially regulated, divided in 446 up-regulated genes (red) and 112 down-regulated genes (green). Genes not 
exceeding a fold change of 1.5 are depicted as basal level (black) at this time point. Results are divided in 5 clusters representing up/down-regulation 
and time of involvement including corresponding GO-BP terms and KEGG pathways. 
doi:1 0.1 371 /journal.pone.01 04548.g002 



tions of IL-6 were detected at three difTerent time points: 2 hours 
(not significant), 4 days, and 7 days p.i. (Figure 7A). Remarkably, 
IL-6 concentrations returned to basal level 6 hours p.i. Production 
of IL-la by epithelial cells was most likely a result of LPS or 
tracheal cytotoxin (TOT) exposure [59]. The enhanced concen- 
tration of IL-6 2 hours p.i. was most likely caused by lung-resident 
immune cells, such as macrophages and mast cells. These cells are 



the firstto recognize bacterial infection [13]. Elevated concentra- 
tions of chemokines CXCL9 (MIG) and CXCLIO (IP- 10) were 
measured between 14-28 days p.i., which suggests a Thl response 
[60]. CCLll (Eotaxin) was enhanced at 14 days p.i. (Figxire 7E). 
Eotaxin is involved in allergic inflammation and is chemoattractive 
for eosinophils [61] and CCRS"^ Th2 lymphocytes [62]. IL-9 was 
increased 7 days p.i. (Figure 7F). IL-9 is mainly produced by T- 
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Figure 3. Principal component analysis of pulmonary gene 
expression. An evolving gene expression profile was indicated after 6. 
pertussis infection of naive mice (IVIean of n = 3), illustrated as principal- 
component analysis (PCA). PCA is a mathematical algorithm [117], 
which describes data based on (dis)similarity. Therefore, a greater 
distance between points in the plot corresponds to a greater 
dissimilarity. In this figure, the similarity of the 10 time points are 
compared based on the expression profiles of the 558 differentially 
expressed genes. Results indicate that 81% and 9% of the variance 
between time points could be addressed by principal component 1 
(PCI) and PC2, respectively. The next principle components contribute 
for the residual variances. The plot revealed a notable change in gene 
expression at 2 days p.i. compared to 4 hours p.i., and largest 
differences at 14 days p.i. The gene expression at 28 days p.i. was still 
far from the naive state. 
doi:1 0.1 371/journal.pone.01 04548.g003 

cells and originally associated with Th2 responses [63]. However, 
there is evidence that IL-9 is produced by a specific T-cell subset: 
Th9 cells [64] and Thl 7 cells [65]. A Thl 7 response was identified 
by markers, such as IL-17A, CCL20 (MIP-3a) and granulocyte 
colony-stimulating factor (G-CSF) [34]. Enhanced concentrations 
of IL-17A and CCL20 were determined in mouse sera 14 days p.i. 
(Figure 7G-H). CCL20 binds exclusively to CCR6, which is 
highly expressed on ThI7 cells [32,33]. G-CSF was produced 
between 4 and 28 days p.i. G-CSF is important for recruitment, 
difierentiation and proliferation of neutrophils [66] and is linked to 
a Thl7 response [67]. No significant changes in concentrations 
were detected for the other twenty-four cytokines. 

In conclusion, serum cytokine profiles indicate the induction of 
an acute phase response (IL-la and IL-6) followed by a multi- 
flavored immune response with Thl (CXCL9 and CXCLIO), Th2 
(CCLl 1 and IL-9) and Thl 7 signatures (G-CSF, CCL20, IL-9 and 
IL-17A). 

Cellular composition of the spleen 

The cellular composition of the spleen was examined (Figure 8A 
and Figure S2) by flowcytometry with a panel of fluorescence- 
labeled monoclonal antibodies binding to ceU-speciflc membrane 
markers (Figure S3 and S4). Results indicate that in the naive 
situation, the majority of splenocytes consist of B-ceUs and T-cells 
(approximately 40% each). Upon B. pertussis infection, the ratio 
between B-ceUs and T-cells changed, since an increased percent- 
age of T-ceUs and decreased percentage of B-ceUs was observed on 



day 2-4 p.i. However, this ratio was restored on day 7 p.i. The 
distribution of minor cell populations in the spleen, such as DCs, 
neutrophils and CD 14"*" cells, changed over time upon B. pertussis 
infection. Here, we observed an early increased percentage of 
CD14"" cells (day 1 p.i.) and DCs (day 1 and day 7-14 p.i.). 
Furthermore, a gradual increased percentage of neutrophils from 
day 4 until day 2 1 p.i. was observed. Importantiy, enlargement of 
the spleen was observed during the course of infection (not 
quantified), suggesting influx or proliferation of ceUs. This might 
explain the observed changes in cellular distribution. 

In conclusion, cellular composition of the spleen indicates 
fluctuations over time upon B. pertussis infection. 

Gene expression in the spleen 

Gene expression in the spleen of infected mice was determined 
by microarray analysis. The study identified 798 genes of the 
whole genome that were diflferentiaUy regulated (|!)-value£0.001, 
FRSl.5): 342 genes were up-regulated, 425 genes were down- 
regulated, and 31 genes that were initially down-regulated and 
subsequently up-regulated (Figure 8B). All involved genes are 
listed in Table S2. A limited set of genes showed altered gene 
expression between 4 hours and 7 days p.i. However, large 
changes in gene expression were observed 21 days p.i., while the 
gene expression was almost back to the original values 28 days p.i. 
Only eight genes still showed altered expression. Over-represen- 
tation analysis (ORA) of the 798 differentially expressed genes in 
the spleen revealed 66 GO-BP terms and two KEGG pathways. A 
selection of 16 terms and 2 pathways is presented in Figure 9, 
which are involved in transcription, ceU cycle, chromatin assembly 
and a number of immunological processes. The corresponding 
proteins are involved in proliferation, difierentiation or activation 
of immune cells in the spleen. Down-regulated genes encode 
proteins that are involved in transcription, adaptive immune 
responses, antigen processing, antigen presentation and chemokine 
signaling (Figure 8B and 10). Down-regulation of these genes 
suggests suppression of the immune response or migration of 
immune cells. 

Immune response. Forty-one immune response genes were 
identified (Figure 10), 13 of which are associated with the adaptive 
immune response. Interestingly, down-regulation of genes in- 
volved in T-ceU activation, B-cell activation, antigen processing 
and presentation, and the subsequent BCR signaling pathway 
were observed 21 days p.i. (Figure 10). Up-regulation of genes 
encoding immunoglobulin heavy chains Ighg, Igha and Igh-vj558 
indicate production of immunoglobulins IgG and IgA by B-cells 
(Figure 10). Down-regulation of the Ighm gene suggests a 
suppressed IgM production or a class-switch to IgG or IgA. 

Membrane receptors. Up-regulation of genes encoding 
myeloid ceU membrane markers, Fc-receptors and membrane 
attack complex (MAC) inhibitory proteins (Cd59a/Cd59h) 
suggests increased presence of myeloid cells, such as neutrophils, 
between 10 and 21 days p.i. (Figure 10). The indurtion of the 
MAC -inhibitory genes suggests inhibition of MAC formation, 
which is part of the complement cascade. In addition, BioGPS 
analysis (Figure S5) revealed the presence of myeloid cell markers 
that are expressed on neutrophils {Trem3 and Clec4e) or 
macrophages (Clec4e). Down-regulated gene expression of mem- 
brane markers of B-cells [Cd22, Cd72, Cd74 and Siglecg) and T- 
cells (Sema4d, Tcra) implies decreased presence of lymphocytes 21 
days p.i. (Figure 10). 

Transcription. Ninety-five genes that were differentially 
regulated in the spleen turned out to be involved in transcription 
(Figure 9), of which 70 genes were down-regulated (data not 
shown). 
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Figure 4. Functional annotation of pulmonary gene expression data after B. pertussis infection. Over-representation analysis (Benjamini- 
corrected p-value£0.05) in GO-BP and KEGG databases of pulmonary gene expression data resulted in respectively 122 terms and 20 pathways 
significantly enriched. A selection of terms and pathways, based on immunological relevance and exclusion of overlapping terms, is given in this 
figure. GO-BP terms were divided in mother terms (dark blue) and corresponding daughter terms (light blue). For each term the amount of genes, the 
percentage of the total term population and the Benjamini-corrected p-value are shown. 
doi:1 0.1 371 /journal.pone.01 04548.g004 



Cell cycle. Most of the genes involved in the cell cycle process 
(31 out of 45) were up-regulated (Figure 10) which suggest 
increased cell replication between 10 and 21 days p.i. in the spleen. 

Regulation of apoptosis. Apoptosis was characterized by 
up-regulation of 12 genes (between 7 and 21 days p.i.) and down 
regulation of 24 genes (21 days p.i.) (Figure 10). The role of 



apoptosis is especially interesting because of the recently described 
link between apoptotic markers and influenza vaccine responsive- 
ness [68]. However, no overlap was found between the markers 
from the study of Furman et al. and the current study. 

In summary, most prominent differences in gene expression in 
the spleen were found around 2 1 days p.i. These differences point 
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Figure 5. Gene expression profiles of pulmonary cytokines and membrane receptors after B. pertussis infection. Results for cytokines 
and membrane receptors include expected source and known corresponding receptors, target cells and receptor ligands. IVIembrane receptors are 
divided in eight different groups according to function. (Mean of n = 3). 
doi:l 0.1 371/journal.pone.01 04548.g005 



to down-regulation of immunological processes. After 28 days p.i. 
the expression profiles returned to basal level, at the same moment 
that the bacteria were almost cleared from the lungs. Interestingly, 
when combining cellular composition (Figure 8A) with gene 
expression profiles (Figure 8B and Figure S5), the increase in 
neutrophils (10-21 days p.i.) was observed with both techniques. 
However, changes in the distribution of B-cells and T-ceUs on day 
2 and 4 p.i. were not detected on gene expression level. Overall, 
these findings suggest that gene expression profiles cannot directly 
be translated in cellular fluctuations in the spleen. 

Overlap of transcriptomic responses in lungs and spleen 

A comparison of transcriptomic datasets was made to determine 
the related biological processes that occurred in both spleen and 
lungs (Figure S6). Seventy-two genes were differentially expressed 
in spleen and lungs, which were divided in four groups: (i) 20 genes 
that were up-regulated in both tissues, (ii) 42 genes that were up- 
regulated in lungs and down-regulated in the spleen, (iii) 2 genes 
that were down-regulated in lungs and up-regulated in the spleen, 
and (iv) 8 genes that were down-regulated in both tissues. Fifteen 
genes of the second group were involved in the immune response, 
such as leukocyte activation, BCR signaling pathway, antigen 
processing and presentation, and chemotaxis. Remarkably, the 
immunoglobuHn-mediated immune response was the only process 
up-regulated in both tissues, suggesting local and systemic 
antibody production. 

In summary, both tissues showed activation of antibody- 
mediated immunity upon a B. pertussis infection. In addition, 
most differentially expressed genes representing immunological 
processes were observed up-regulated in the lungs and down- 
regulated in the spleen. 

Pertussis-specific IglVl, IgG and IgA antibody responses 

Antibody isotypes, IgM, IgG and IgA against B. pertussis were 
analyzed for a period of 28 days p.i. The formation of IgM against 
B. pertussis antigens present in outer membrane vesicles (OMV) 
was not observed during the course of infection (Figure 1 lA). 
Contrarily, anti-B. pertussis IgG in mouse sera was detectable 
from 10 days p.i. and the responses were progressively increased 
until 28 days p.i. (Figure IIB). Notably, no IgG titers were 
detectable against the individual purified B. pertussis antigens 
pertactin (Prn), pertussis toxin (Ptx), filamentous hemagglutinin 
(FHA) or fimbria type 2 and 3 (Fim2/3) (Figure S7). However, 
anti-OMV IgG was significantly increased at 14, 21 and 28 days 
p.i. (Figure IIC) and consisted of IgGl, IgG2a and IgG2b 
subclasses, but not IgG3. (Figure IID-G), Remarkably, IgGl 
levels dropped significantly at 28 days p.i. compared to 2 1 days p.i. 
In addition, significantly increased anti-OMV IgA responses were 
measured in sera (Figure 1 IH) as well as in lung lysates 
(Figure 1 II) 28 days p.i. 

In conclusion, intranasal infection leads to a broad spectrum of 
systemic antibody subclasses against outer membrane proteins of 
B. perlnssis within 14 clays p.i., but responses against purified 
virulent antigens Ptx, Prn, FHA and Fim2/3 remain undetectable 
until 28 days p.i. Moreover, B. pertussis -speciSic anti-OMV IgA 
antibodies were secreted in the lungs. 



CD4^ T-cell responses after B. pertussis infection 

Cytokine profiles of antigen-specific memory 004"*" T-ceUs in 
the spleen were determined 66 days p.i. Splenocytes of naive and 
infected mice were stimulated for 8 days with Ptx, FHA or Prn, 
important virulence factors present in the acellular pertussis 
vaccine. Cytokine profiling on single cell level was established by 
using intracellular cytokine straining (ICS), while cytokine 
concentrations in splenic culture supernatants were determined 
by using a MIA. CD44 expression was used to focus on activated 
CD4^ T-cells. Results of ICS indicate a significant increase in 
IFNy- and IL- 1 7 A-producing Prn-, FHA-, and Ptx-specific 
CD4"^CD44'' T-cells after infection (Figure 12A and B). Except 
for FHA-stimulation, no evidence was found for specific IL-5- 
producing CD4 CD44 T-cells (Figure 12C). Cytokine profiles in 
culture supernatant revealed a significant increase in TNFa 
production by splenocytes of infected mice after FHA stimulation 
and IL-17A production after Prn, Ptx, or FHA stimulation 
(Figure 12D-F). No changes were observed in IFNy production in 
culture supernatants (Figure 12D-F). Notably, FHA stimulation 
induced besides Thl and Thl7 cytokines an increase in IL-5, IL- 
13, IL-4, and IL-10 production by splenocytes of infected mice 
(Figure 12E). IL-5 and IL-13 production occurred also after Prn 
stimulation (Figure 12D). This might indicate that these three B. 
pertussis antigens induce a mixed Th-subset CD4"'" T-ceU response, 
of which Thl and Thl 7 are the major components. While ICS 
specifically focuses on CD4''"CD44''" T-ceUs, multiple cell types can 
produce the accumulated cytokines present in the splenic culture 
supernatants after stimulation. The production of IL-17A is 
exclusive for Thl 7 cells [69]. Besides Thl cells, IFNy could also be 
produced by CTLs and NK-cells. The high levels of IFNy 
produced by splenocytes of naive mice after Ptx stimulation may 
be explained by TLR4 activation on innate cells by Ptx [70]. A 
high background of Prn-responding IL-17A CD4 CD44 T-cells 
and IL-17A in supernatant was seen in two out of six naive mice, 
indicating Pm-activated IL-17A production in unprimed spleno- 
cytes. 

In summary, CD4 memory T-cell analysis revealed that 

infection with B. pertussis primarily generates specific Thl and 
Thl 7 responses, regarded essential in protection against a B. 
pertussis infection [12,20,21]. 

Discussion 

An integrated systems biology approach was applied to 
investigate the evolving immune response of mice following a 
primary intranasal B. pertussis infection, considered capable of 
inducing sterilizing immunity. Data sets from individual analytical 
platforms, such as microarray, flow cytometry, multiplex immu- 
noassays, and colony counting, addressing different biological 
samples taken at multiple time points during the immune response 
were merged. Altogether, a comprehensive overview of immuno- 
logical processes occurring upon primary infection of mice with B. 
pertussis was constructed in three consecutive phases, as shown in 
Figure 13: (i) the innate phase, (ii) the bridging phase, and (iii) the 
adaptive phase. 

Highlights featured in this study are: (i) a protracted acute phase 
response, as part of the innate phase, lasting for at least 28 days, 
which co-occurred with prolonged antigen exposure; (ii) the 
recognition of several B. pertussis 'PAMPs' (pathogen associated 



PLOS ONE I www.plosone.org 



8 



August 2014 I Volume 9 | Issue 8 | e104548 



Molecular Signatures of a Bordetella pertussis Infection 



o 

CM T- 

Q Q Q Q Q 




LCN2 
SLPI 
CHI3L1 
CXCL17 




Antimicrobial peptides 



0.05 0.13 0.33 0.671.5 



.> o -"l- T- oo 

Z^TQQQQQQQ 



o 00 
CM ■>J T- CM (N 

Q Q Q Q D Q Q 




^1 



Acute phase 




RETNLA 
SAA3 
LRG1 
TNFAIP2 

STEAP4/TNFAIP9 

ITIH4 

TFF2 

MMP3 

TIMP1 

SERPINA3N 

SERPINA3G 

SERPINA3F 

LTF 

ADM 

SAA4 

0RM1 

PTX3 

GM15056/DEFB52 
A2M 



CTSK 
IRAK2 
ZBP1 / DAI 
IRF7 

CLEC10A/MGL1 
CD209F/SIGNR8 
CD209G / SIGNR7 
CLEC4A3 / DCL3 
IRF5 

LY86 / MD-1 




Phagocytosis 




Antigen presentation 
MHC-I 



o 00 
CM 1- t- CM CM 

D D D D Q Q Q 



Pathogen recognition 



D 



TFDP2 
NFIL3 

SFPil /PU.1 

IRF7 

STAT1 

JAK3 

IRF8 

IRF5 

SYK 

PTPN6 / SHP-1 



Transcription factors 




CIITA 

PSMB8 

TAP1 

GD74 

H2-AA 

H2-AB1 

H2-EB1 

H2-DMA 

H2-DMB1 

H2-DMB2 

FGER1G 

FCGR2B 

UNC93B1 

IFI30 

CD4 

HSPA1A 

HSPA5 



Antigen presentation 
MHC-II 




C3AR1 



Complement 
Classical pathway 



O •<- CO 

CM 1- CSi CM 

Q Q Q Q Q Q Q 




G3AR1 



Complement 
Alternative pathway 



CM 

Q Q 



o ^ 

h- CM 

Q Q Q Q 



I LILRB3 
FCGR2B 
IFITM1 
BLNK 
RAC2 
PIK3GD 
PLCG2 
CD22 
CD79A 
PTPN6 
VAV1 
CARD11 
GRB2 
PRKCB 
CD72 



BCR signaling 



o 

CM T- t- CM 

Q Q Q Q Q Q 




Immunoglobulins 



Figure 6. Pulmonary gene expression representing immunological pathways as function of time after B. pertussis infection. Details on 
pulmonary gene expression profiles indicated expression of anti-microbial peptides, acute phase proteins, pathogen recognition, phagocytosis, both 
classical and alternative pathways of complement cascade, transcription factors, antigen presentation via both IVIHC-I and MHC-II, B-cell receptor 
(BCR) signaling and immunoglobulin related genes. (Mean of n = 3). 
doi:1 0.1 371 /journal.pone.01 04548.g006 
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Figure 7. Serum cytokine concentrations as function of time after B. pertussis infection. The concentrations of (A) IL-lot, (B) lL-6, (C) CCL1 1, 
(D) CXCL9, (E) IL-17A, (F) IL-9, (G) CXCL10, (H) CCL20, and (I) G-CSF in serum were analyzed over time. Results indicated that after the acute phase 
response (IL-1 -x and IL-6) a multi-flavored immune response with Thi (CXCL9 and CXCL1 0), Th2 (CCL1 1 and IL-9) and Thi 7 (G-CSF, CCL20, IL-9 and IL- 
17A) was induced. Data represented as mean concentrations including individual values (n = 3). p-values were determined by one-way ANOVA with 
multiple comparison compared to naive mice (* = p<0.05, ** = p<0.01, *** = p<0.001 and **** = p<0.0001). 
doi:1 0.1 371/journal.pone.01 04548.g007 



molecular pattern), i.e. LPS, DNA and GalNac, (iii) early 
signatures for recruitment of adaptive immune cells e.g. the 
cytokine CCL20 that attracts CCR6^Thl7 cells, (iv) expression of 
specific membrane markers (e.g. DAP12) and transcription factors 
that determine the strength and direction of the adaptive immune 
response. Furthermore, features of adaptive immunity consisted of 
(v) IgA secretion in the lungs; (vi) a broad systemic antibody 
response; (vii) a Thi and Thl7-mediated response in the lungs and 
(viii) systemic formation of specific memory Thi and Thi 7 cells. 

Various signatures, associated with the induced protective 
immune response as presented in Figure 13, are discussed in 
detail hereafter. The innate phase was characterized by the 



presence of AMPs, acute phase proteins and chemo tactic cytokines 
as detected by microarray analysis and MIA (Figure 5—7). 
Expression of the involved genes is most likely initiated by 
activation of TLR2 and TLR4 by LPS of fi. pertussis [8,10]. Up- 
regulated gene expression of four AMPs {CM311 , Cxcll 7, Lcn2, 
Spli) was found in the lungs during the whole course of infection. 
These AMPs contribute to innate protection as was shown for 
other respiratory pathogens [38-40,42]. In addition, these AMPs 
have the capacity to recruit immune cells, such as neutrophils, 
DCs, and monocytes [37,41] and are therefore important for the 
induction of subsequent immune responses [36]. 
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Figure 8. Splenic cellular composition and gene expression profiles as function of time after B. pertussis infection. (A) The distribution 
of the splenic cellular composition illustrated as pie charts at a selection of time points during the 6. pertussis infection. FACS was used to determine 
the percentage (mean of n = 3) of B-cells (CD19*), CD4* T-cells, CDS* T-cells, neutrophils (Gr-1*), macrophages (GDI lb*), dendritic cells (33D1*), NK- 
cells (DX5*) and CD14* cells (CD14*) in splenocytes. (B) Fold changes of expression were calculated compared to naive mice and significant gene 
expression results (FR>1.5, P-value<0.001) were visualized as a heatmap (Mean of n = 3). This test was performed for each gene at each time point. In 
total, 798 genes were differentially regulated (FRal.5, p-value<0.001) compared to naive mice, divided in 342 up-regulated (red) and 425 down- 
regulated genes (green). Additionally, a distinct group of 31 genes showed down-regulation between 4 hours and 4 days p.i., but were up-regulated 
between 10 and 21 days p.i. Genes not exceeding a fold change of 1.5 are depicted as basal level (black) at this time point. Results are divided in 3 
clusters representing up/down-regulation and time of involvement and included corresponding GO-BP terms and KEGG pathways. 
doi:10.1371/journal.pone.0104548.g008 



Furthermore, the expre.ssion of particular cytokines on gene and 
protein level in the innate phase (Figure 5 and 7), in combination 
with increased gene expression of cytokine receptors and the 
change in cellular composition of the spleen (Figure 5 and 8A), 
gave detailed insight into the immunological processes during the 
course of infection, as depicted in Figure 13. For example, the 
gene expression levels of the cytokines Cxcll, Cxcl2, Cxcl5 and 
Cxcll 7 was in agreement with the recruitment of neutrophils and 
DCs towards the site of infection [41,71,72]. Furthermore, the 
gene expression of the cytokines Cxcl9, Cxcll 0, Cxcll 1 and Ccl20 
was increased early after infection. CXCL9, CXCLIO and 
CXCLl 1 are known to be involved in the recruitment of 
CXCR3* Thl cells, and CCL20 in the recruitment of CCRe"" 
Thl7 cells [32,33,60]. The up-regulation of Ccr6 and Cxcr3 genes 
was observed at day 14 p.i. in the lung (Figure 13). 



Remarkably, gene expression of cytokines, acute phase proteins 
and AMPs was temporary suppressed at 2 days p.i. Simultaneous- 
ly, the proliferation of bacteria started in the lungs. These results 
indicate temporary suppression of the immune system by B. 
pertussis. Several immune evasion strategies of B. pertussis have 
been described in literature [73-77]. The evasion is attributed to 
virulence factors, such as adenylate cyclase toxin (ACT), pertussis 
toxin (Ptx), tracheal cytotoxin (TCT) and filamentous hemagglu- 
tinin (FHA). These investigators demonstrated enhanced produc- 
tion of IL-10 and regulatory T-ceUs, delay of neutrophil 
recruitment, and inhibition of early chemokine production in 
mice in relation to secretion of 5. pertussis virulence factors. In the 
present study, these immunosuppressive properties of 5. pertussis 
might explain the suppression of the differential gene expression 
observed 2 days p.i. In addition, secretion of OMVs occurs during 
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Figure 9. Functional annotation of splenic gene expression data after B. pertussis infection. Over-representation analysis (Benjamini- 
corrected p-value<0.05) of all 798 differentially expressed genes in the spleen in GO-BP and KEGG databases, This resulted in respectively 66 terms 
and 2 pathways significantly enriched of which a selection, based on immunological relevance and exclusion of overlapping terms, of 1 6 terms and 2 
pathways are presented. GO-BP terms were divided in mother terms (dark blue) and corresponding daughter terms (light blue). For each term the 
amount of genes, the percentage of the total term population and the Benjamini-corrected p-value are shown. 
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Figure 1 0. Splenic gene expression profiles of specific processes as function of time after B. pertussis infection. Expression profiles are 
summarized for immune response, cell cycle, regulation of apoptosis, T-cell activation, B-cell activation, immunoglobulin, B-cell receptor (BCR) 
signaling, antigen processing, cytokines and membrane receptors in spleen. (Mean of n = 3). 
doi:l 0.1 371/journal.pone.01 04548.g01 0 



bacterial colonization [78]. These OMVs are thought to serve as 
decoy for the host immune response and their immunomodulatory 
properties support survival of the bacteria [79,80] . 

During the bridging phase, the number of bacteria in the lungs 
increased (Figure 1). The highest number of bacteria in the lungs 
was found at 7 days p.i. followed by a decrease-, suggesting that day 
7 is the turning point in the progress of the infection. 
Simultaneously, microarray analysis suggests the presence of 
neutrophils, macrophages, and DCs in the lungs based on BioGPS 
analysis (Figure SI). These cells are likely responsible for (i) 
phagocytosis, (ii) pathogen recognition, (iii) antigen presentation 



and (iv) differentiation by expression of cell receptors in the period 
between 4—10 days p.i. (Figure 13). 

The innate immune cells appeared to cause the initial decrease 
in numbers of bacteria in the lungs from 7 days p.i., because a 
strong antibody response (IgM, IgG, IgA) (Figure 1 1 ) and CD4''' 
cell response were not yet present. Based on gene expression of 
neutrophil markers, the neutrophil invasion in the lungs probably 
occurred from 4 days p.i., which is observed on ceUular level by 
McGuirk et al. [81]. This influx of cells is most likely orchestrated 
by IL-17A and G-CSF [66,67,69]. In the present study, both 
cytokines were enhanced in the sera of infected mice. Further- 
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Figure 11. Antibody profiling after B. pertussis intect'ion. Antibody responses were determined In mouse sera after intranasal Infection by 
whole cell 6. pertussis ELISA or MIA and expressed as OD450 or Fluorescent Intensity (F.I.) respectively. (A) IgM against OMV B1917 were absent. (B) 
Whole cell 6. pertussis ELISA Indicated IgG antibody formation 10-14 days p.i., which Increased until day 28 days p.i. (C-G) Levels of total IgG and 
Individual subtypes (IgGI, lgG2a, lgG2b and lgG3) against OMV B1917 Indicated presence of multiple subclasses. (H) Antl-OMV IgA antibodies In 
serum were Induced after 14-28 days p.i. (I) In lung lysates, antl-OMV IgA antibodies were detected 14 and 28 days p.i. I p-values were determined by 
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more, genes involved in phagocytosis were found to be up- 
regulated (Figure 6) which is important for protection against B. 
pertussis [55,82]. 

This study demonstrated also that B. perlussis infection activates 
different pattern recognition receptor (PRR) pathways, such as 
TLR4, TLR9, DAI and several CLRs (Figure 6). Bacterial ligands, 
e.g. LPS, DNA and GalNAc trigger these signaling pathways 
[8,45,47]. Up-regulation of genes involved in TLR4-signalling 
[Md-l, IrfS) is probably caused by detection of bacterial LPS 
[8,10,43]. Bacterial DNA recognition occurred through the TLR9 
pathway [Clsk, Irf7) and the cytosolic DNA sensor (Dai). The 
receptor DAI recognizes both host DNA, upon tissue damage, and 
bacterial DNA, which leads to an innate immune response [45]. 
Furthermore, increased gene expression of the CLRs, Signr7, 
SignrS and Mgll was observed, which are predominandy 
expressed on DCs [28,83]. The enhanced expression of mouse 
DC-SIGN-related genes Signr? and SignrS during the B. perlussis 
infection is a novel finding. Despite indications of sensing of 6- 
sulfo-sialyl Lewis(x) oUgosaccharide by Signr 7, the exact ligands 
and role of these receptors remains to be discovered [46]. 
Recognition of GalNAc via Mgll has recently been associated 
with mast cell activation by -6. pertussis [47]. An important 
conclusion is that ligands of 5. pertussis activate multiple pathogen 
recognition pathways during the infection. This may be important 
for the induction of a broad adaptive immune response, including 
Thl as well as Thl7 cells, mucosal IgA and systemic IgG. 

During the bridging phase, up-regulation of MHC class I and II 
genes were measured, enabling antigen presentation to both CD4"'" 
Th-cells and cytotoxic CD8^ T-cells (CTLs). Several studies 
focused on the role of Th-c:ells because- of their important 
contribution to protect against B. pertussis [12,20,21]. Neverthe- 
less, the CTLs seem also important for protection against B. 
pertussis, indicated by CTL infiltration in lungs and formation of 
FHA-specific IFNa-producing CTLs upon B. pertussis infection as 
was previously found in mice and humans [81,84]. 

In addition, the bridging phase was characterized by simulta- 
neous gene expression of severed receptors on innate immune cells 
in the lungs and spleen. Interestingly, this included MAIRs, 
TREMs, semaphorins and DAP12, which is the first evidence that 
these receptors might play a role in B. pertussis infection. 
Especially DAP12, which is abundantiy expressed in lungs, 
interacts with other membrane receptors and influences cytokine 
expression [85]. Coupled to MAIR-II, DAP12 is an inhibitor of B- 
ceU responses [86]. Together with TREM2, DAP12 down- 
regulates TLR and FcR expression [87]. Overall, these former 
receptors function as an important bridge between innate and 
adaptive immune responses, determining the strength and 
direction of the adaptive immune response by the successive 
activation of lymphocytes, such as CTLs, Th-cells, and B-ceUs 
[29,88-91]. Since the exact role of these receptors in the 
protection against B. pertussis remains unresolved, additional 
research is needed to find whether targeting of these receptors 
might be a viable strategy to elicit desired immune responses. In 
summary, the bridging phase suggests the presence of APCs and 
myeloid cells, such as DCs and neutrophils based on gene 
expression of membrane markers. However, this was not vahdated 
by histological data. 

Finally, the adaptive phase was recognized by the formation of 
T-cell and B-cell responses and a broad humoral response 
(Figure 13). In this same period, the naive mice were able to 
clear B. pertussis from their lungs. Gene expression, cytokine 
secretion and cellular analysis revealed presence of activated CD4 
Th-cells. The highest number of CD4^ Th-ceUs was detected at 14 
days p.i., which is in agreement with a study showing the largest 



influx of Thl 7 cells in the lungs after B. pertussis infection at 14 
days p.i. [12]. Phenotyping of CD4'^ Th-cells, based on cytokine 
profiles in serum (protein) and lung (genes) and expression of 
specific transcription factors or membrane markers, suggests a 
broad Th-response in the present study. The Thl response in the 
lungs was characterized by Ifng and Stall expression [92] . Gene 
expression of the mucosal homing receptor Ccr6, Cardll, Cell 9, 
Ccl21a, Ccl21c, and signature cytokines III 7a and III 7f might 
indicate Thl 7 cells in the lungs [32,69,93-96]. Upon re- 
stimulation of splenocytes with purified B. pertussis antigens, 
Prn- and FHA-specific Thl and Thl 7 t^pe CD4^ T-cells could be 
detected in the memory phase (66 days p.i.) (Figure 12). In 
addition, enhanced IL-9 levels in serum suggest activation of Th9 
cells and Thl 7 cells as these cells produce large amounts of IL-9 
[64,65]. However, Tregs, NKT-ceUs and mast cells are also 
sources of IL-9. In summary, the different phases of T-cell 
development were mapped; the expansion phase of naive T-cells 
in the spleen (2 days p.i.), the effector phase in the lungs (14 days 
p.i.), contraction phase in the spleen (2 1 days p.i.) and finally the 
specific memory T-cells (66 days p.i.) (Figure 5, 8A, 10 and 12). 
Data from this study suggests generation of local specific Thl/ 
Thl 7 cells in the lungs. Thl cells, which produce IFNy are 
important for clearantx of S. pertussis. Postponed clearance was 
shown in IFNa-depleted mice [97]. Thl 7 cells have been indicated 
as key cells to control bacterial infections at mucosal sites [98]. 
Recent studies revealed that Thl 7 cells increase the clearance of 
B. pertussis after an intranasal infection in animals [12,13,99]. 

Infection with B. pertussis induces also B-cell and aiitiliody 
mediated immunity. Microarray analysis showed increased gene 
expression of the BCR in the lungs. In addition, gene expression of 
IgM, IgG and IgA was observed in lungs and spleen. Whereas, 
most genes for antibody isotypes were found up-regulated in both 
tissues, IgM was down-regulated in the spleen (Figure 6 and 10). 
This was confirmed by the detection of IgA and IgG (IgGl, IgG2a 
and IgG2b) after 14 days p.i. However, IgM antibodies were not 
found systemicaUy (Figure 1 1). In mice, IgGl is associated with a 
Th2-like response, while IgG2a (in Balb/c mice) and IgG2b 
suggest induction of a Thl response [100]. IgG2b is also linked to 
Thl 7 lymphocytes [101]. AU antibody responses in the present 
study were directed against whole [:ell B. pertussis and outer 
membrane vesicles, but not against the B. pertussis antigens that 
are typically present in acellular vaccines: Ptx, Prn, FHA and 
Fim2/3. However, memory Thl /Thl 7 cells were detected upon 
re-stimulation with Prn and FHA (Figure 12). The polymeric 
immunoglobulin receptor (pIgR) was highly expressed during the 
whole course of fi. pertussis infection in the lungs. pIgR is essential 
for the transport of IgA into the mucus and is seen as a bridge 
between innate and adaptive mucosal responses [102,103]. The 
expression of pIgR was drastically enhanced 14 days p.i., leading 
to secretion of mucosal IgA in the lungs (Figure 5 and 1 1). 
Interestingly, the increased pIgR expression coincided with IL- 
1 7A production in the lungs and sera (Figure 5 and 7), which is in 
agreement with the finding that Thl7-mediated responses 
influence the local humoral response by inducing pIgR expression 
and elevating secretory IgA levels [104]. Therefore, the transport 
of IgA to the mucosa, which is orchestrated by IL-17A via the 
induction of pIgR, supports an important role for local immunity. 
In summary, B. pertussis infection induces a broad humoral 
response, which was observed by gene and protein expression. 

Our study also contains several conceptual points for discussion. 
First, altered gene expression in the spleen caused by B. pertussis 
infection was modest, except for 21 days p.i. Two possible reasons 
can be given: (i) the infection initiated a local immune response in 
the lungs and (ii) gene expression is required to be rather strong to 
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Figure 12. Systemic memory Th-cytokine profiles. (A-C) Splenocytes of naive and infected mice were collected 66 days p.i. and in vitro 
restimulated with Prn, FHA or Ptx for 8 days. The percentages of IFNy-, IL-17A-, and lL-5-producing CD3*CD4*CD44* T-cells were determined using 
ICS. (D-F) Cytokine levels in supernatant after 7 days of stimulation were determined by using a MIA. Results (mean of n = 6) are corrected for the 
background level in the presence of medium as control. Statistical differences between the groups were detected for the ICS with a non-parametric 
IVIann-Whitney test and for the MIA with a Student t-test on the log-transformed data. ** = p<0.01, *** = p<0.001, **** = p<0.0001. 
doi:1 0.1 371 /journal.pone.01 04548.g01 2 



exceed the significance threshold. After the lungs and the draining 
lymph nodes, the spleen is the third organ in line involved in the 
generation of the immune response. Therefore, gene expression 
effects could be less abundant. Furthermore, the spleen consists of 
many cell types with different gene expression profiles, so changes 
in expression levels in a single cell type might not have been 
detected. Cell sorting of individual cell types and subsequent gene 
expression profiling could overcome this problem. Second, specific 
antibody and T-cell responses were only measured for a limited set 



of available purified antigens. The inclusion of other antigens, such 
as outer membrane proteins, would perhaps have allowed a more 
detailed analysis. Third, in this study we propose murine infection- 
induced immune signatures as a benchmark to be used in the 
development of improved human pertussis vaccines. Since B. 
pertussis is not a natural pathogen for mice, translation of the 
results obtained in this study to the human situation remains to be 
interpreted with caution. Nevertheless, using a (naive) murine 
model has several advantages. It is regarded an important small 
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Figure 13. Overview of immune response evolving after B. pertussis infection. Merging data from a wide array of functional analyses reveals 
a comprehensive overview of the immune response of mice upon 6. pertussis infection, (i) Innate phase, (ii) bridging phase and (ill) adaptive phase are 
indicated in the response facing variable levels of 6. pertussis present in the initial inoculum, the exponential growth phase, and in the clearance 
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animal model for pertussis vaccine development [105] and enables 
to dissect the local and systemic immune response in great detail, 
as shown in this study, rather than using blood or nasal washes 
from humans undergoing a B. pertussis infection. Furthermore, 
the role of individual gene signatures and products could be 
assessed in future research in mice by using knockin or knockout 
strains, or by functional interventions. Encouraging is that many 
observations made in this study, such as the rapid clearance upon 
reinfection and involvement of Th 1 7 and Th 1 responses, are in 
line with recent studies on B. pertussis infection in baboons, which 
resemble the human situation more closely [13,15]. Translation of 
biomarkers from mouse to man however should be performed 
with care. 

Overall, our study describes the molecular and cellular sequence 
of events that eventually lead to the induction of protective 
immunity to pertussis, including local immunity in the lungs, 
which is not induced by current vaccines. Especially the inductions 
of Th 1 and Th 1 7 cells and IgA in the lungs are regarded key 
elements in superior immunity. Our data can be used to select 
vaccine concepts that resemble infection, with respect to e.g. 
administration route and prolonged antigen exposure, and 



immune signatures elicited. Local delivery of vaccines may be an 
important advancement for improved protection as was previously 
demonstrated with a live attenuated pertussis vaccine [106]. Also 
adjuvants could be selected and optimized that can steer host 
immunity of inactivated vaccines towards an infection-like 
immune response with minimal adverse effects. The molecular 
and cellular fingerprint of the local immune response discovered in 
B. pertussis infected mice thus can provide a strong guidance for 
developing improved pertussis vaccines. 

Methods 

Ethics Statement 

The independent ethical committee for animal experimenta- 
tions 'Dierexperiment Commissie (DEC)' of the National Institute 
for Public Health and the Environment (RIVM) reviewed the 
animal experiments in this study according to the guidelines 
provided by the Dutch Animal Protection Act. The committee 
approved documents with the identification numbers 
'DPA20 1100230' and 'DPA201 100348'. 
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Materials 

Bordetella pertussis strain B1917 [107] was isolated in The 
Netherlands in 2000 and was kindly provided by Frits Mooi. 
Bordet Gengou agar plates with 15% sheep blood were purchased 
from BD (Gat. No. 254400, BD, The Netherlands). Verweij 
medium was obtained from Bilthoven Biologicals (BBio, Bilthoven, 
The Netherlands). PBS (pH 7.2) and Roswell Park Memorial 
Institute Media 1640 (RPMI) were purchased from Invitrogen 
(Gibco, Invitrogen). All RPMI media used in this experiment were 
supplemented with I0"/o fetal calf serum (FGS, Hyclone, New 
Zealand), 100 units penicillin, 100 units streptomycin, and 
2.92 mg/ml L-glutamine (p/s/g; Invitrogen). Shockbuffer for 
splenocyte isolation consisted of 8.3 g/L NH4CL, 1 g/LNaHCOs 
and 5000 lE/L Heparine dissolved in dHgO (pH 7.4; ice cold) and 
was filtered (0.22 |xm) before use. FAGS buffer consisted of PBS 
pH 7.2 supplemented with 0.5% BSA (A3803, Sigma Aldrich, 
Germany) and 0.5 mM EDTA (IGN Biomedicals). Horseradish 
peroxidase (HRP) labelled goat-anti-mouse IgG and i?-Phycoer- 
ythrin (RPE)-conjugated goat-anti-mouse IgM, IgA, IgG, IgGl, 
IgG2a, IgG2b and IgG3 were obtained from Southern Biotech 
(Southern Biotech, United States) Fluorescently labeled antibodies 
for flow cytometry, anti-GD3-Pacific blue (17A2), anti-GD4- 
Pacific blue (GK1.5), anti-IA-IE-PerCP-Cy5.5 (M5/1 14.15.2), 
anti-33Dl-APG (33D1) and anti-IL-5-APC were purchased from 
Biolegend (ITK, The Netherlands). Fluorescendy labeled antibod- 
ies, anti-GD4-APG (RM4-5), anti-GD8a-PerGP (53 6.7), anti- 
GD69-FITG (HI. 2 F3), anti-Gr-l-PE (RB6-8G5), anti-GD19- 
FITG (1D3), anti-CD14-PE (rmG5-3), anti-DX5-APG (DX5), 
anti-GDI Ib-PE (Ml/70), anti-GD40-FITG (HM40-3), anti- 
CD44-FITG (IM7) and anti-IFNy-PE were purchased from BD 
(BD Biosciences, The Netherland.s). Fluorescently labeled anti- 
bodies, GDllc - PE-Gy7 (N418) and anti-IL-17a-PerGP-Gy5.5 
were obtained from eBioscience (eBioscience, Austria). Live/dead 
staining was purchased from Invitrogen (LIVE/DEAD Fixable 
Aqua Dead Cell Stain Kit; Invitrogen). Anti-CD 16/CD32 was 
purchased from BD (BD, The Netherlands). 

Cultivation of Bordetella pertussis 

Bordetella pertussis strain B1917 was grown on Bordet-Gengou 
agar plates for 4 days at 35 °C. Subsequently, multiple colonies 
were transferred onto new plates and grown for 24 hours at 35°G, 
which was repeated 3 times. Colonies on each plate were 
harvested in 1.5 ml Verweij medium [108] and the suspensions 
were pooled. After centrifugation (30 min., 1932 g), cells were 
washed with 15 ml Verweij medium. Subsequently, bacterial 
suspensions of approximately 2xlO'" cfu/ml were resuspended in 
Verweij medium containing 15% glycerol (v/v) and aliquots of 
0.25 ml were flash frozen in ethanol and dry ice and stored at — 
80°C. 

Challenge culture 

A stock of B. pertussis B1917 was diluted with Verweij medium 
to a final concentration of 5x10' cfu/ml. The number of cfu was 
confirmed by plating 100 |iL of the infection suspension (1:2000 in 
Verweij medium) on Bordet-Gengou agar plates. Plates were 
incubated for 4 days at 35°C and cfu were counted by using a 
colony counter (ProtoCOL, Synbiosis, Cambridge, United King- 
dom). 

Animal experiment 

Female BALB/ c mice (Harlan, The Netherlands), 8-weeks-old, 
were divided in groups of three animals and housed in cages 
(macrolon III including filter top). Three healthy mice were 



euthanized at day 0 and considered as naive group. Other mice 
were intranasally infected under anesthesia (isoflurane/oxygen), 
with 2x10^ cfu B. pertussis B191 7 in 40 |a,l Verweij medium at day 

0. Groups of animals (n = 3) were euthanized after 2, 4, 6 hours, 
and after 1, 2, 4, 7, 10, 14, 21 and 28 days of infection, 
respectively. At 56 days p.i. five groups of primary infected animals 
and five groups of naive animals (n — 3) were infected as previously 
described and euthanized after 4 hours, and after 2, 7, 10 and 14 
days of infection. In addition, six infected and six healthy mice 
were euthanized at day 2 1 and day 66 for the analysis of GD4^ T- 
cell r(;sponsf;s. Micx' w(;re bkxi, under anesthesia (isoflurane/ 
oxygen), by orbital bleeding and sacrificed by cervical dislocation. 

Whole blood from each mouse was collected in a blood 
collection tube (MiniGoUect 0.8 ml Z Serum Sep GOLD, Greiner 
Bio-One, Austria). After cxiagulation (10 min. at room tempera- 
ture), sera were taken after centrifugation (10 min., 3000 g) and 
stored at — 80°C. Lung lobes were excised. The right lobe was 
placed in 900 (il Verweij medium and kept at room temperature. 
The left lobe was placed in 1 ml RNAlater (Qiagen), incubated 
overnight at 4°C and stored at — 80°C. The spleen was excised 
and divided in two equal parts. One pic-cc was placed in 5 ml of 
RPMI medium and kept on ice. The other piece was placed in 
1 ml RNAlater (Qiagen), incubated overnight at 4°C and stored at 
— 80°C. For the analysis of 004""" T-ceU responses, whole spleen 
was collected and placed in 5 ml of RPMI medium and kept on 
ice. 

Colonization assay 

Lung tissue in 900 jtl Verweij medium was homogenized for 
15 seconds using a Labgen 7 Homogenizer (Cole Palmer, 
Schiedam, The Netherlands). Homogenizer was sanitiz(;d bc-tween 
samples with ethanol (80%) and water. Lung homogenates were 
serially diluted (undiluted, 1:10, 1:100 or 1:1000) in Verweij 
medium depending on the expected number of cfu. Suspensions of 
100 |xl were plated on Bordet-Gengou agar plates and incubated 
for 4 days at 35°C. The numbers of colonies were counted by 
using a colony counter (ProtoCOL, Synbiosis, Cambridge, United 
Kingdom) and calculated as cfu per mouse. 

RNA isolation 

Lung and spleen tissue was homogenized by using a Labgen 7 
Homogenizer (Cole Palmer, Schiedam, The Netherlands) in 700 
and 2100 |i,l Qiazol (Qiagen Benelux, Venlo, The Netherlands), 
respectively. RNA isolation was performed using a miRNeasy 
Mini Kit with DNAse treatment (Qiagen Benelux, Venlo, The 
Netherlands) according to the manufacturer's protocol. RNA 
concentrations were determined by UV spectroscopy (Tech3 
module. Synergy Mx, BioTek, Winooski, United States). RNA 
quality was determined by using electrophoresis (RNA nano 6000 
kit, 2100 Bioanalyzer, Agilent Technologies, Amstelveen, The 
Netherlands). Results were given as RNA integrity numbers (RIN) 
between 1 nd 10 (manual Bioanalyzer). Samples with a RIN of 7 
or higher were used for microarray analysis. For microarray 
analysis of lung tissue, RNA concentrates of individual mice were 
analyzed for the following time points: naive, 4 hours, 2 days, 4 
days, 7 days, 10 days, 14 days, 21 days and 28 days post infection 
(p.i.). For spleen tissue, samples of individual mice were analyzed 
for the following time points: naive, 4 hours, 1 day, 2 days, 4 days, 
7 days, 10 days, 14 days, 21 days and 28 days p.i. 

Microarray analysis 

Amplification, labeling and hybridization of RNA samples to 
microarray chips was carried out at the Microarray Department of 
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the University of Amsterdam, The Netherlands, as described 
previously [109]. 

Briefly, 500 ng total RNA of each sample was amplified (aRNA) 
according to the Agilent QuickAmp kit manual (Agilent technol- 
ogies). A common reference sample for either lung or spleen tissue 
was made by pooling equimolar amounts of aRNA from 
individual samples of the respective tissues. Cy3 and Cy5 mono- 
reactive dyes (GE Healthcare) were used to label individual 
samples and the common reference sample, respectively. Labeling 
was performed with 10 ml of CyDye solution and incubated for 
1 hour before the reaction was quenched by adding 5 ml 4 M 
hydroxylamine (Sigma- Aldrich). Purification was performed by 
using a clean-up kit (E.Z.N.A. MicroElute RNA Clean Up Kit, 
Omega Bio-Tek, Norcross, United States). The yields of amplified 
RNA and incorporation of CyDye were determined by using a 
spectrophotometer (NanoDrop ND-1000, NanoDrop products, 
Wilmington, United States). 

Hybridization mbcture was prepared by adding a dried mixture 
(1:1) of sample (Cy3) and common reference (Cy5) together with 
sample tracking control (STC, Roche NimbleGen) and hybrid- 
ization cocktail (NimbleGen Arrays User's Guide -Gene Expres- 
sion Arrays Version 5.0, Roche NimbleGen). Samples were 
ini ul)atcd successively for 5 min. at 95°C and 5 min. at 42°C. 
Hybridization was performed by loading a sample onto a 
microarray (NimbleGen 12x135 k Mus musculus, Roche, Ger- 
many) containing probes for 44,170 genes with 3 spots per target 
probe. Hybridization was performed with a NimbleGen Hybrid- 
ization System 4 (Roche NimbleGen) for 20 hours at 42°C. After 
washing (NimbleGen Arrays User's Guide - Gene Expression 
Arrays Version 5.0), slides were scanned in an ozone-free room 
with a microarray scanner (Agilent DNA microarray scanner 
G2565CA, Agilent Technologies). Each microarray corresponded 
to labeled RNA from specific tissue (lung or spleen) from one 
individual mouse. Feature extraction was performed with Nim- 
bleScan \'2.5 (Roche NimbleGen) resulting in a table containing 
individual probe signal intensities for both dyes. Complete raw and 
normalized microarray data and their MIAME compliant 
metadata from this publication have been submitted to the 
GEO database (www.ncbi.nlm.nih.gov/ geo) and assigned the 
identifier GSE53294. 

Data analysis of gene expression 

Quality control was performed on raw data by means of Cy3- 
Cy5 scatter plots, and by comparing signal average and 
distribution across slides. All slides passed quality control. Raw 
microarray data for gene-coding probes were normalized in R 
(www.r-project.org), by using a four step approach [109]: (1) 
natural log-transformation, (2) quantile normalization of all scans, 
(3) correcting the sample spot signal for the corresponding 
reference spot signal and (4) averaging data from replicate probe 
spots. Further analysis of normahzed data of 44170 probes was 
performed in R and Microsoft Excel. 

Genes difiFerentially expressed between experimental groups 
(naive and various time points p.i.) were identified by using 
ANOVA. Fold ratio induction or repression of individual genes 
was calculated by comparing mean gene expression levels of 
infected groups to the group with naive mi[:e. Data are presented 
as the average normalized gene expression levels of three mice per 
group. Probes were considered differentially expressed if they met 
the following two criteria: (i) a p-value<0.001 (ANOVA), which 
corresponds to a Benjamini-Hochberg False discovery rate (FDR) 
[110] of <5%; and (ii) an absolute fold ratio >1.5 (infected 
compared to naive mice) for at least one time point. If multiple 
probes corresponding to the same gene were significant, their data 



were averaged to remove redundancy for further analysis. Gene 
expression of each group was compared with the group of naive 
mice. Differences in gene expression were visualized in heat maps, 
(GeneMaths XT, Applied Maths, St-Martens-Latem, Belgium). 
Genes were arranged according to similar expression patterns in 
time at which genes exceeded the fold ratio cut-off of 1.5. To 
facilitate visual interpretation of heat maps, only induction (red) 
and repression (green) of gene expression levels with fold ratios > 
1.5 are visualized, therefore presenting fold ratios <1.5 as naive 
level (black). 

Functional enrichment with an over-representation analysis 
(ORA) was carried out by using DAVID [1 1 1] based on Gene 
Ontology' Biological Processes (GO-BP) and Kyoto Encyclopedia 
of Genes and Genomes (KEGG). Cell-specific or tissue-specific 
gene sets were obtained by using data from BioGPS [112,113]. For 
each cell type or tissue, gene expression data were compared to the 
average of aU tissues included in the data set. The top 100 genes, 
which had the highest tissue-specific expression and were included 
on the NimbleGen array, were selected. If the BioGPS dataset 
included multiple entries for the same cell type (e.g., 004"*" and 
CD8^ T-ceUs), the resulting sets were combined in our subsequent 
analysis. 

Wliole cell B. pertussis ELISA 

Immunoglobulin G responses against B. pertussis were deter- 
mined in sera by an ELISA. Therefore, B. pertussis B1917 was 
heat-inactivated for 45 min. at 56°C, after which ELISA plates 
(Immunolon 2HB, Thermo Scientific) were coated overnight at 
room temperature with 100 |jl/well whole cell B. pertussis 
(OD59o„m = 0.1 in PBS). Plates were washed with 0.03% (v/v) 
Tween 80 in water. Then, plates were incubated for 1 hour at 
37°C with 100 nl/well serial dilutions of serum samples in PBS, 
pH 7.2 (Gibco) with 0.1% (v/v) Tween 80. Plates were washed 
twice and incubated (1 hour, 37°C) with 100 |J,l/well of Goat anti 
Mouse IgG-HRP (1:5000) in PBS with 0.1% (v/v) Tween 80 and 
0.5% (w/v) Protifar (Nutricia). After washing twice, 100 nl/weU 
peroxidase substrate (4.2 mM Tetra Methyl Benzidine and 
0.012% H2O2 in 0.11 M sodium acetate buffer, pH 5.5.) was 
added, plates were incubated at room temperature and the 
rea[:tion terminated after 10 min. I)y adding 100 |Xl/well, 2 M 
H2SO4 (BBio, Bilthoven, The Netherlands). Finally, the absor- 
bance was recorded at 450 nm with a plate reader (BioTek reader 
EL808, Bio-Tek, USA) and antibodies against whole cell B. 
pertussis are presented in OD450 signal, because a reference is not 
available. 

Multiplex immunoassay (MIA) for antibody response 

Responses of IgM, IgA, total IgG and the IgG subclasses (IgGl, 
IgG2a, IgG2b and IgG3) against B. pertussis antigens P.69 
pertactin (P.69 Prn), filamentous hemagglutinin (FHA), pertussis 
toxin (Ptx), combined fimbria type 2 and 3 antigens (Fim2/3) and 
outer membrane vesicles B1917 (OMV B1917) were determined 
in sera by using a MIA. Outer membrane vesicles from B. pertussis 
B1917 were produced as previously described [114] with 
additional changes. Conjugation of antigens and OMVs to beads 
was performed as described previously [115]. Serum was diluted 
1/100 in PBS containing 0.1% Tween 20 and ?>"/a bovine serum 
albumin and mixed 1:1 with conjugated beads. After incubation 
with i?-Phycoerythrin (RPE)-conjugated anti-mouse IgM (1:100), 
IgA (1:100), IgG (1:200), IgGl (1:200), IgG2a (1:40), IgG2b (1:200) 
and IgG3 (1:200), samples were analyzed by using a Bio-Plex 
system (Bio-Plex 200, BioRad). Total IgG titers (U/ml) against B. 
pertussis antigens were calculated from mean fluorescent intensity 
(MFI) by using a reference from the National Institute for 
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Biological Standards and Control (Code 99/520, NIBSC, 
England) in the Bio-Plex Manager software (Bio-Rad Laborato- 
ries). Results for IgM, IgG, IgG subclasses and IgA antibodies 
against OMV B1917 are presented in fluorescent intensity (F.I.), 
because a reference is not available. 

Lung lysate preparation and IgA antibody assay 

The lung lysates used for the colonization assay were 
centrifuged (2000 g, 10 min.) and filt(;r(;d (0.22 \im Millex-GV, 
Merck KGaA, Darmstadt, Germany) to remove cell debris and 
bacteria. Filtrates were stored at — 80°G. Prior to analysis, lung 
lysates were diluted 1:10 in PBS containing 0.1% Tween 20 and 
3% bovine serum albumin. Analysis of IgA antibodies in the lung 
lysate was performed against B. pertussis antigens (PtK, Prn, FHA, 
and Fim2/3) and OMV B1917 as described in the section 
multiplex immunoassay (MIA) for antibody response. 

Splenocyte isolation 

Spleen tissue was forced through a 70 |im cell strainer (BD 
Falcon, BD Biosciences, USA) and washed two times with 5 ml 
medium (RPMI, FCS, p/s/g) to obtain a single cell suspension. 
After centrifugation (550 g, 5 min.), 1 ml shock buffer was added 
to cell pellet to lyse erythrocytes. After one minute, shock 
incubation was stopped by adding 25 ml medium. For cell count 
(CASEY model TT, Innovatis, Roche, Germany), cells were 
centrifuged and resuspended in 1 ml medium. 

Cellular composition of spleen 

Splenocytes (1x10^) were transfered to each well of a V-bottom 
96-wells plate. After centrifugation (400 g, 3 min.), pellets were 
washed with 100 (tl FACS buffer. Fc-specific receptors on cells 
were blocked with anti-CD 16/CD32 in FACS buffer (2%, v/v) by 
5 min. incubation on ice to prevent aspecific binding of labeled 
antibodies to cells. Cells were stained after centrifugation by using 
a mixture (100 |tl) of labeled antibodies (30 min., 4°C). To 
determine the number and identity of cells, three diSerent panels 
of labeled antibodies were used to prevent overlapping fluores- 
cence spectra. Panels consisted of 1) total T-cell population (CDS), 
CD4 T-cell (CD4), CD8 T-cell (CDS), early activated T-cell 
(CD69), and neutrophil (Gr-1); 2) T-cell (CDS), B-cell (CD 19), 
monocyte (CD 14), and natural killer cell (DX5); 3) T-cell (CD3), 
macrophage (CDl lb), DC (33D1, CDl Ic), MHC class II (lA/IE), 
and activation marker (CD40). Data acquisition was performed 
with a FACSCanto II (BD Biosciences, USA) and data analysis 
was done by using Flo\vJo software (Tree Star Inc., USA). Debris 
and dead cells were excluded by using forward/ side scatter 
characteristics and live/ dead staining, respectively. 

Stimulation of splenocytes 

Splenocytes were resuspended in the culture medium IMDM 
(Gibco) +8% FCS+p/s/g+20 |a,M P-mercaptoethanol (Sigma) and 
cultured (24-wells plate; 6x10'' cells/well) for 7 days at 37 °C in a 
humidified atmosphere containing 5% CO2. Culturing was in the 
presence of either the culture medium, or the culture medium 
supplemented with Prn P. 69 (in house expressed and purified as 
described previously [1 16]), heat-inactivated (15 min at 95''C) Ptx 
(Kaketsuken, Japan) or FHA (Kaketsuken, Japan) at a final 
concentration of 1 |ig/ml. Supernatants were collected after 7 
days for cytokine analysis. Subsequently, cells were transferred to 
U-bottom 96-wells plates (5x10'^ cells/well) and stimulated 
overnight by using the same antigens and conditions and used 
for intracellular cytokine staining. 



Intracellular cytokine staining 

The intracellular cytokine staining was performed on 8-days 
stimulated splenocytes by using a Cytofix/Cytoperm Fixation/ 
PermeabHization Solution Kit (BD Biosciences, USA). Briefly, cells 
were incubated during the last 5 hours of stimulation with 10 Hg/ 
ml Golgiplug (BD Biosciences), 1 ^g/ml aCD28 (BD Pharmin- 
gen), and 1 Jig/ml aCD49d (BD Pharmingen). Cells were stained 
after washing with FACS buffer with anti-CD4-Pacific blue, anti- 
CD44-FITC and live/ dead staining. Subsequently, the cells were 
fixed, permeabUized, and stained with anti-IFNy-PE, anti-IL-5- 
APC and anti-IL-17a-PerCP-Cy5.5. Data from fluorescence- 
activated cells were acquired on FACS Canto II (BD Biosciences) 
and analyzed with Flowjo software (Tree Star Inc., USA). CD44 
was used to distinguish between naive and antigen experienced T- 
cells. IL-5 was used as a Th2 cytokine, IFNy as a Thl cytokine, 
and IL-17 as a Thl 7 cytokine. 

Multiplex immunoassay (MIA) for cytokines 

Cytokine (IL-4, IL-5, IL-10, IL-13, IL-17A, TNFa, and IFNy) 
concentrations (pg/ml) were determined in culture supernatants of 
7 days stimulated splenocytes of individual mice by using a MIA 
(Milliplex mouse cytokine 7-plex luminex kit; Merck KGaA, 
Darmstadt, Germany). Cytokine (Eotaxin, G-CSF, GM-CSF, 
IFN-y, IL-10, IL-12 (p40), IL-12 (p70), IL-13, IL-15, IL-17A, IL- 
la, IL-ip, IL-2, IL-3, IL-4, IL-5, IL-6, IL-7, IL-9, IP-10, KC, LIE, 
LIX, M-CSF, MCP-1, MIG, MlP-la, MIP-ip, MIP-2, RANTES, 
TNF-a, VEGE) concentrations (pg/ml) present in serum were 
determined by using a MIA (MiUiplex MAP Mouse Cytokine/ 
Chemokine - Premixed 32 Plex; Merck KGaA, Darmstadt, 
Germany) according to the manufacturer's instructions. Briefly, 
serum/supernatant (25 |a,l) was diluted in supphed assay buffer (1/ 
1; v/v), mixed with 25 (il conjugated beads and incubated 
overnight at 4°C on a plate shaker. Subsequently, beads were 
washed twice with washing buffer by using a Bio-Plex handheld 
magnetic washer (BioRad, USA) and incubated with 25 (il 
detection antibodies (1 hour, 20-25 C) followed by incubation 
with 25 H-l streptavidin-phycoerythrin (30 min., 20-25°C). After 
washing twice with 200 jtl washing buffer, beads were resuspended 
in 150 ill PBS. Samples were analyzed using a multiplex system 
(Bio-Plex 200, BioRad, USA). In total, 50 beads per region were 
analyzed in 100 |il sample volume. Mean fluorescence intensity 
(MFI) of each cytokine was converted into pg/ml by using a 
dilution series of a standard (3.2-10,000 pg/ml) and Bio-Plex 
Manager software 5.0 (Bio-Rad, USA). 

MIP-3a (CCL20) cytokine ELISA 

Serum concentration of Macrophage Inflammatory Protein 3a 
(MIP3a/CCL20) was determined by ELISA (Mouse MIP-3a 
ELISA kit, RAB006 1 , Sigma- Aldrich, Germany) according to the 
manufacturer's manual. In short, serum (50 |J,1) was diluted in 
assay buflTcr A (1/1; v/v) and incubated in an antibody-coated 
ELISA plate (overnight, 4°C, gentie shaking). After washing four 
times with 300 |J,1 wash buffer, wells were incubated with 100 |a.l 
biotinylated antibody (1 hour, room temperature, gentle shaking). 
After washing four times with 300 |jl wash buffer, weUs were 
incubated with 100 |J,1 streptavidin solution (45 min., room 
temperature, gentie shaking). After washing four times, weUs were 
incubated with 100 ml TMB one-step substrate reagent (30 min., 
RT) and reaction stopped with 50 ml stop solution. FinaUy, the 
absorbance was recorded at 450 nm with a plate reader (BioTek 
reader EL808, Bio-Tek) and converted to MIP-3a concentration 
by using a dilution series of a standard (2.06-1,500 pg/ml) and 
Gentech 5 software (BioTek). 
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Statistical analysis 

Cellular composition of spleen, antibody titers and cytokines 
levels were analyzed by using a one-way ANOVA with multiple 
comparisons of all time points post infection vs. naive mice 
followed by a Dunnett t-test. ICS data was analyzed by using a 
Mann- Whitney t-test. For cytokine profiling after re-stimulation of 
splenocytes, MIA data was first log-transformed before t-test 
analysis. In all corresponding figures, j!<-values are represented as 
*=^<0.05, **=p <0.01, ***=^<0.00I and ****=^<0.0001. 

Supporting Information 

Figure SI Cell type comparison analysis for lung gene 
expression. Data represent gene profiles in the lung per cell type 
or tissue extracted from BioGPS. Results from B-cells, mast cells, 
T-cells, dendritic cells, neutrophils and lung are depicted. (Mean of 
n = 3). 
(TIF) 

Figure S2 Cellular composition of the spleen as function 
of time after B. pertussis infection. The percentage of (A) 
CD 14"" cells (CD 14"*"), (B) NK-cells (DX5+), (C) macrophages 
(CDl lb"*), (D) neutrophils (Gr-1"^, (E) T-ceUs (CD3"", CD4'"CD8" 
and CD8"'CD4"), (F) B-ceUs (CD 19"*) and (G) dendritic ceUs 
(33D1"*) in splenocytes were analyzed over time, ^-values were 
determined by one-way ANOVA with multiple comparison 
compared to naive mice (*=^<0.05, **=jb<0.01, ***=p< 
0.001 and ****=jf)<0.0001). (mean of n = 3). An increased 
percentage- of CD 14"*" cells (monocytes, DCs and macrophages) 
was found 1 day p.i., followed by natural killer (NK) cells and 
macrophages at 2 days p.i. There was a gradual increase in the 
percentage of neutrophils in the splci-n imtil 7 days p.i. 
Furthermore, an increased percentage of T-c(-lls, distinguished 
by the CDS marker, was found 2 and 4 days p.i. A similar increase 
was observed for both CD4+CD8^ and CD8*CD4" T-cell 
subsets, indicating that this increase was not specific for either 
T-helper or cytotoxic T-ceUs. The CD4"*" cells increased signifi- 
candy 21 days p.i. and remained constant until at least 28 days p.i. 
B-cells decreased during two periods: 2-4 and 10-14 days p.i. 
Finally, a significant increase in percentage of DCs (33D1"*) was 
detected 7 days p.i. 
(TIF) 

Figure S3 Flow cytometric-gating strategy for CD 14+, 
Gr-1+, CD19+ and CDllb+ cell population analysis in 
splenocytes. Each FACS plot shows results from one mouse of 
the group. Per cell type, the change in numbers of cells is visible 
over time for CD 14+, Gr-1+, CD 19+ and CDllb+ cell 
populations. Time points 2 hours, 6 hours and 10 days p.i. were 
excluded in this figure. (*) FACS plot for CD 14+ c(;lls at da)' 4 and 
day 7 p.i. were manually adapted by BiExponential transforma- 
tion in Flowjo for better visualization without influencing the 
gating. 
(TIF) 

Figure S4 Flow cytometric-gating strategy for CD3+, 
CD4+, CD8+, DX5+ and 33D1+ cell population analysis 
in splenocytes. Each FACS plot shows results from one mouse 
of the group. Per cell type, the change in numbers of cells is visible 
over time for CD3+, CD4+, CD8+, DX5+ and 33D1+ cell 
populations. Time points 2 hours, 6 hours and 10 days p.i. were 
excluded in this figure. (*) FACS plot for 33D1+ cells at day 1, day 
7 and day 14 p.i. were manually adapted by BiExponential 
transformation in Flowjo for better visualization without influ- 
encing the gating. 
(TIF) 



Figure S5 Splenic gene expression profiles of cell types 
and tissues extracted from BioGPS databases. Data 
represent gene profiles of the following cell t^pes or tissues: B- 
cells, neutrophils, mast cells, T-cells, macrophages, dendritic cells 
(DCs), lung, lymph nodes and spleen. (Mean of n = 3). 
(TIF) 

Figure S6 Genes differentially regulated in both lung 
and spleen. Comparison of transcriptomic data in lung and 
spleen revealed an overlap of 72 genes differentially regulated in 
both tissues. Genes divided in 4 groups based on expression 
pattern. Functional annotation showed that genes in group 2 were 
mostly involved in immunological processes. (Mean of n = 3). 
(TIF) 

Figure S7 IgG antibody response in serum against Ptx, 
Prn, FHA and Fim2/3 after B. pertussis infection. IgG 

antibody titers against four purified B. pertussis antigens; (A) 
pertussis toxin (Ptx), (B) pertactin (Prn), (C) filamentous hemag- 
glutinin (FHA), and (D) Fimbriae 2 and 3 (Fim2/3) were 
determined using a multiplex immunoassay (mean of n = 3). The 
antibody titers were determined in mouse sera 14, 21 and 28 days 
after an intranasal infection with B. pertussis. No antibodies were 
detected for Ptx, Prn and Fim2/3. For FHA, 1 out of 3 mice 
showed induced titers antibody titers at 21 days p.i. 
(TIF) 

Table SI Full list of genes expressed in the lungs. This 
file contains the fuU list of significantly expressed genes obtained by 
microarray analysis of the lungs, which can be used for additional 
in-depth information for each gene. In total 558 genes were 
significandy altered compared to naive mice. The expression 
pattern for each gene is given including the fold changes at each 
time point. In addition to the official gene symbol (obtained from 
the Nimblegen array), the full gene name, synonyms and gene ID 
are given. Synonyms are not given for all genes but only for genes 
that were included in detail in the manuscript. Gene ID is a linked 
to gene information on the NCBI website. Subsequently, there are 
three columns with gene information obtained in databases; (i) 
BioGPS, (ii) Gene Ontology - Biological Processes (GO-BP), (iii) 
KEGG pathways. Pull-down menus can be used to search for 
specific information in the fuU gene list. For BioGPS, information 
about cell and tissue specific gene expression is given for the 
following databases; B-cells (B), dendritic cells (DC), lung (Lu), 
lymph node (LN), mast cells (MC), macrophages (MF), neutrophils 
(N), spleen (Sp) and T-ceUs (T). For GO-BP the following 
biological processes were included; 'acute inflammatory response', 
'acute-phase response', 'antigen processing and presentation', 'cell 
activation', 'chemotaxis', 'defense response', 'immune response', 
'inflammatory response', 'intracellular signaling cascade', 'phago- 
cytosis', 'regulation of apoptosis', 'response to bacterium' and 
'transcription'. For KEGG the following pathways were included; 
'Antigen processing and presentation', 'Apoptosis', 'B ceU receptor 
signaling pathway', 'Cell adhesion molecules (CAMs)', 'Chemo- 
kine signaling pathway', 'Complement and coagulation cascades', 
'Cytokine-cytokine receptor interaction', 'Endocytosis', 'Leukocyte 
transendothelial migration', 'MAPK signaling pathway', 'NOD- 
Uke receptor signahng pathway' and 'ToU-like receptor signaling 
pathway'. 
fXLSX) 

Table S2 Full list of genes expressed in the spleen. This 
file contains the fuU list of significantly expressed genes obtained by 
microarray analysis of the spleen, which can be used for additional 
in-depth information for each gene. In total 798 genes were 
significandy altered compared to naive mice. The expression 
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pattern for each gene is given including the fold changes at each 
time point. In addition to the official gene symbol (obtained from 
the Nimblegen array), the full gene name, synonyms and gene ID 
are given. Synonyms are not given for all genes but only for genes 
that were included in detail in the manuscript. Gene ID is a linked 
to gene information on the NCBI website. Subserju(;ntly, there are 
three columns with gene information obtained in databases; (i) 
BioGPS, (ii) Gene Ontology - Biological Processes (GO-BP), (iii) 
KEGG pathways. PuU-down menus can be used to search for 
specific information in the full gene list. For BioGPS, information 
about cell and tissue specific gene expression is given for the 
following databases; B-cells (B), dendritic cells (DC), lung (Lu), 
lymph node (LN), mast cells (MC), macrophages (MF), neutrophils 
(N), spleen (Sp) and T-ceUs (T). For GO-BP the foUowing 
biological processes were included; 'acute inflammatory response', 
'antigen processing and presentation', 'cell activation', 'chemotax- 
is', 'chromatin assembly', 'defense response', 'hemopoiesis', 
'immune response', 'inflammatory response', 'intrac:ellular signal- 
ing cascade', 'regulation of apoptosis', 'response to bacterium' and 
'transcription'. For KEGG the following pathways were included; 
'Antigen processing and presentation', 'Apoptosis', 'B cell receptor 
signaling pathway', 'Cell adhesion molecules (CAMs)', 'Chemo- 
kine signaling pathway', 'Complement and coagulation cascades'. 
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